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ABSTRACT 

Early-type galaxies obey a narrow relation traced by their stellar content be- 
tween the mass and size (Mass- Radius relation). The wealth of recently acquired ob- 
servational data essentially confirms the classical relations found by Burstein, Ben- 
der, Faber, and Nolthenius, i.e. logi?!/ 2 oc logM^ 054 for high mass galaxies and 
logi?!/ 2 oc logM^ ' 3 for dwarf systems (shallower slope), where R\/2 and M s are the 
half-light radius and total mass in stars, respectively. Why do galaxies follow these 
characteristic trends? What can they tell us about the process of galaxy formation? 
We investigate the mechanisms which concur to shape the Mass-Radius relation, in 
order to cast light on the physical origin of its slope, its tightness, and its zero point. 
We perform a theoretical analysis, and couple it with the results of numerical hy- 
drodynamical (NB-TSPH) simulations of galaxy formation, and with a simulation of 
the Mass-Radius plane itself. We propose a novel interpretation of the Mass- Radius 
relation, which we claim to be the result of two complementary mechanisms: on one 
hand, the result of local physical processes, which fixes the ratio between masses and 
radii of individual objects; on the other hand, the action of cosmological global, sta- 
tistical principles, which shape the distribution of objects in the plane. We reproduce 
the Mass-Radius relation with a simple numerical technique based on this view. 
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1 INTRODUCTION 

Galaxies are known to exist in a variety of morphological 
types going from dwarf objects (in turn grouped in irregu- 
lars, ellipticals and spheroidals), to disk galaxies with the 
central bulges of different size (the origin of which is not yet 
clearly established), and finally to the elliptical ones of large 
mass and dimensions. The mass and size among the different 
types may vary by orders of magnitudes. This taxonomy im- 
plies a coherent picture of galaxy formation and evolution in 
cosmological context. Indeed, the situation observed in the 
local Universe is further complicated going back into the 
past at higher and higher redshifts, where a large variety of 
objects seems to exist. 

Galaxy formation in cosmological context. In the 
presently accepted model of the Universe, dominated by cold 
Dark Matter (DM) and by a mysterious form of energy ex- 
pressed by the cosmological constant, and containing Bary- 
onic Matter (BM) and photons (thereinafter A-CDM Uni- 
verse), the cosmic structures are originated from the grav- 
itational collapse of DM halos that gives rise to complexes 
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of larger and larger scale within which baryons (gas) infall, 
forming stars and galaxies in a complicate game of star for- 
mation (SF), chemical enrichment, gas heating and cooling, 
and galactic winds. Many years of observations and theoreti- 
cal speculation have clarified that building up galaxies of dif- 
ferent morphological type (schematically irregulars, spirals 
and ellipticals) or even different parts of the same galaxy, 
e.g. the bulge, halo and disc in spirals, cannot be reduced 
to a unique physical mechanism. While there is sort general 
consensus on the formation of spiral and irregular galaxies, 
the question is unsettled in the case of early-type galaxies 
(ETG), the massive ones in particular. It is unclear how the 
A-CDM cosmological model should be reconciled with the 
observational claims for large, a nd red, galaxies alrea dy in 
place at very high redshifts (e.g. lHarrison et al.ll201ll ). The 
problem is twofold: on one hand, it must be explained how 
and when star formation is quenched in both massive and 
small haloes, which is necessary to reconcile the theoreti- 
cal prediction with the observed galaxy mass function (see 
e.g. iBundv et al1l2006f ); on the other hand, one should also 
clarify how such massive systems as those recently observed 
can form at very high redshifts. Historically, two opposite 
scenarios exist: 
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(1) the monolithic view, which predicts that most of the 
mass in massive spheroids is already assembled in a strong 
and rapid burst of SF at high rcdshift. Galaxies evolve al- 
most passively thereafter, without being strongly affected 
by su bsequent mergers (see[Ch 10S1 fe Carraro 2002; P eebles! 
120021 . and references). This picture does not conform to the 
present theoretical, concordance scenario of cosmology; 

(2) the hierarchical view, according to which DM haloes 
(and BM inside) grow by mergers of smaller units in the 
dry (gas-poor, no SF ) or wet (gas-rich and SF) mode (see 
|Pe Lucia et al.l 120061 . and references therein). However the 
role played by mergers (major, minor, wet and dry) is still 
uncertain (|Hopkins et al.l 120101 ) and several different ver- 
sions of the same tell are around (e.g. iGonzalez et al.ll20lj ). 
Detailed numerical simulations have shed doubts on the co- 
herence with which mergers can bring galaxies along the 
tight structural re lations observed in the local Universe (e.g. 
iNfpoti et ailbOOgD . 

The early hierarchical or quasi monolithic view. The mod- 
ern picture, lately emerging from observational data and 
theoretical investigations, is as follows: at a certain (high) 
redshift, perturbations made of DM and BM detach them- 
selves from the Hubble flow, and collapse on their own form- 
ing a proto-galaxy rich of sub-structures inside the com- 
mon gravitational potential well. They merge, form stars 
and eventually give origin to a single galaxy. Depending 
on the initial density and /or angular momentum, the end 
product is different. Systems with low (or nearly zero) an- 
gular momentum will end up as spherical or elliptical ob- 
jects. The massive ones and/or those originated by strong 
density perturbations had an early dominant episode of SF 
ever since followed by quiescence, whereas the low mass 
and/or lower density ones undergo a series of SF episodes 
along the Hubble time (galactic breathing ). This view is 
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Likely objects with high angular momentum end up as disc 
galaxies, in which the local low density environment favours 
SFs prolonged all over the galaxy life. The bulge (if present) 
likely follows the scheme envisaged for ETGs. Finally, irregu- 
lar galaxies should originate from very low density perturba- 
tions within a shallow potential well undergoing intermittent 
and perhaps delayed SF. The advantage with this scenario 
is that ETGs in place at redshift greate r than 2 are possibl e 
jKobavashill2"005l ; iMerlin fc Chiosill2006l ; iMerlin et al.ll2012h , 
and many structura l and chemica l properties are explained 
(sec lMatteucci|[2007l ; IChiosill2007h . 

Masses and Radii of galaxies. In recent years, much 
attention has been paid to the Mass-Radius Relationship 
(MRR) of galaxies, in particular the ETGs and the compact 
and passive ones at high z. The MRR is indeed basic to 
any theory of galaxy formation and evolution. The subject 
of the MRR of galaxies from ETGs to dwarf ellipticals and 
dwarf spheroidals, including also bulge s and Glob ular Clus- 
ters has been recently reviewed bv lGrahaml |201lf ) to whom 
we refer for many details. The current MRR for ETGs will 
be presented below in great detail. 

In addition to this, convincing evidence has been gath- 
ered that at relatively high redshifts, objects of mass com- 
parable to that of nearby massive galaxies but smaller di- 
mensions exist. These "compact galaxies" are found up to 



z > 3 with stellar masses from 10 to 10 Mq and half-light 
radii from 0.4 to 5 kpc (i.e. 3 to 4 times than more com- 
pact than the local counterparts of the same mass), and in 
nearly similar proportions there are ga laxies with the same 
mass but a variety of dim ensions (e.g.. iMancini et al.ll200^ ; 
Valentinuzzi et all bOlOal ). and bulge to disk ratios (e.g., 



Karim et al.ll2011h . However, we will consider here only the 
case of standard ETGs, leaving the compact galaxies aside 
but for a qualitative suggestion about their interpretation. 

In the present study we ask ourselves the question: Why 
do ETGs obey a rather narrow mass-size relation instead of 
scattering around showing a broader combination of these 
two parameters? Spurred by this, we look for general phys- 
ical principles governing this important scale relation. To 
clarify the aims and the methods of this study, we antici- 
pate here the essence of our analysis. We speculate that the 
observed MRR for ETGs is the result of two complemen- 
tary mechanisms. On one hand, the mass function of DM 
haloes hosting the visible galaxies gives (i) the typical cut- 
off mass at which, at any redshift, haloes become "common" 
on a chosen spatial scale, and (ii) the typical epoch at which 
low mass haloes begin to vanish at a rate higher than that 
at which they are born, because of merger events. On the 
other hand, these constraints define two loci (curves) on the 
MR-plane, because to each mass and formation redshift a 
typical dimension (i.e., radius) can be associated (using a 
basic relation between mass and radius of a collapsing ob- 
ject). If the typical dimension of a galaxy is somehow related 
to that of the hosting DM halo (as our NB-TSPH models 
seem to suggest), then the region of the MR-plane between 
the two limits fixed by the halo mass function is populated 
by galaxies whose dimensions are fixed at the epoch of for- 
mation, and only those objects that are "possible" at any 
given epoch may exist, populating a narrow region of the 
MR-plane. 

The paper is subdivided as follows. In Section [5] we 
present and discuss the observational MRR for a sample of 
ETGs galaxies taken from the SDSS catalogue. In Section 
[3] we shortly describe the NB-TSPH models of ETGs that 
we used in the present analysis. In Section [4] we present our 
interpretation of the MRR, based on elementary theories 
of cosmology and galaxy formation. In Section [5] we sim- 
ulate the MRR highlighting the deep causes governing the 
MRR. In Section[5]we advance a possible explanation for the 
anomalous position of compact galaxies at high redshift. Fi- 
nally, in Section [7| we draw some remarks and conclusions. 



2 THE OBSERVATIONAL MASS-RADIUS 
RELATION 

An impressive body of data have been acquired concerning 
the masses and the dimensions of galaxies. T h e situ ation 
has been recently di scussed by IShankar et al.l (|201lh and 
iBernardi et alj (|201lT ) and references therein. 

The normal early-type galaxies. The observational 
data we are conside ring is the HB sample selected by 
IBernardi et"afl (2010) from the SDSS catalogue, containing 
~ 60, 000 galaxiesiT 

The observational MRR is displayed in 
Fig. □ The linear best fit of the SDSS data is 

1 The selection conditions are fracde[v r ]=l and b/a > 0.6, 
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Figure 1. The log[i?]y2] versus log[M s ] relation for galaxies, 
where R1/2 is the half-mass radius, and logM 3 is the total stellar 
mass. Although -R1/2 is not strictly identical to the effective ra- 
dius R e , they are very close to each other. Throughout this paper 
we will always use -R1/2, which is easier to calculate for dynam- 
ical NB-TSPH models of galaxies, and assume R1/2 S Re- The 
contour lines display the regions of the MR-plane populated by 
the same number of objects (the bar yields the correspondence 
between color-cod e and number of galax ies). The data are from 
the HB sample of iBernardi et ah The green solid line is 

the linear best fit of eqn.JT}. 



logifya = 0.54 logM 3 - 5.25 



(1) 



where M s in Mq is the estimated stellar mass (in Mq) and 
R1/2 ( m kpc) is the radius containing half of it (nearly iden- 
tical to the classical effective radius R e ). The slope (and zero 
point) of the above MRR is quite rob ust as it coincid e s with 
previous determinations: using the IBurstein et al.l Jl997 ) 
data, the same slope has been foun d by IChiosi fc Carrard 
|2002l ). and using the SDSS data bv lShen et al l (|2003r i. 

The distribution of th e bulk of galaxi e s is al so confirmed 
by the smaller sample of IShankar et all (|201 lh always ex- 
tracted from the SDSS survey but using slightly different 
selection criteria. The area cove red by the obser v ation al data 
is slightly larger than that the IBernardi et all j2010l ) data, 
fn any case no significant differences can be noted. Our anal- 
ysis of the MRR will stand only on SDSS sample, thus se- 
curing homogeneity of mass and radius estimates. However 
to complete the scene it is worth looking at the position on 
the MR-plane of dwarf galaxies. Since no attempt is made 
to homogenize the data for these latter objects with those 
of the SDSS sample, dwarf galaxies will not be included in 
the analysis of the MRR and will be considered separately. 

The dwarf galaxie s. The dwarf galaxies are taken from 
IBurstein et al.l l| 19971 ). The sample is made of dEs and 



therefore the samp le is dominated by elliptical galaxies (see 
IBernardi et al.ll20ld . for all details). 



dSphs. It is worth recalling that the masses used by 
IBurstein et al.l (|l997l ) are the dynamical masses and not the 
stellar masses, so this group is not strictly homogeneous with 
the sample for ETGs . In addition to the few dwarfs of the 
IBurstein et al.l (| 1997t ) sample, to complete the picture, we 
also consider the mean relationship for the dwarf galaxies 
of the Local Grou p according to the measurements made by 
IWoo et al.l (|2008T l who yield the relationship 



log R 1/2 = 0.3 log M s - 2.7. (2) 

to which also the lBurstein et all l|l997l ) dwarfs seem to obey. 

Changing slope of the MRR. It is soon evident that there 
is no unique slope for the MRR of the different groups of 
objects. The slope is 0.54 for ETGs and 0.3 and lower for 
dwarf galaxies. Looking at the data in detail, the slope is 
even steeper than 0.54 in the region of the largest and most 
massive ETGs go ing up to 1 and e ven m o re, see the top par t 



of the MRR by IBernardi et al.l < |2010h . iGuo et al 
van Dokkum et all l|2010l )," and Fig. 1 of iGrahaml 



2009), 



2011) 



This is a point to keep in mind when interpreting the obser- 
vational data. 

General Remarks. Information and details on how the 
stellar masses M s , and half-mass radii, R1/2, have been 
derived can be found in the original sources to which the 
reader should refer. Of course some possible systematic bi- 
ases among the different sets of data are to be expected, 
whose entity, however, ought to be small. This is some- 
what sustained by the overall agreement among different 
sources as far as some general relationships are concerned, 
e.g. t he agreement in the s lope of the MRR fo r ETG s be- 
tween HemMdreTal] (|20ld ) and lBurstein et all j 19971 ). The 
same for the dwarf galaxies. Perhaps somewhat larger un- 
certainties are present in the case of compact galaxies the 
data of which will be presented in Section [6] However, since 
the three groups of objects will be treated separately and 
only from a general qualitative point of view, no homoge- 
nization of the data is needed. Our discussion is limited to 
the SDSS sample, the data of which is not only internally 
consistent but also homogeneous to those of the theoretical 
simulations. As far as compact galaxies is concerned, what 
will matter here is the general consensus that galaxies of 
very different dimensions and comp arable stellar masses are 
found (e.g. see lShankar et al1l201ll for general discussion of 
the subject). 

Galaxy counts in the MR-plane. Before proceeding further 
it is worth looking at the number frequency distribution of 
galaxies with given mass and radius. We limit ourselves to 
the SDSS sample. To this aim we divide the MR-plane in 
a grid of square cells with dimensions of 0.05 in units of 
A log Ai" s /(10 12 Mq) and Alog-Ri/2, and count the galaxies 
falling in each cell. In Fig. [2]we display the 3D space of these 
parameters to simultaneously highlight the distributions of 
the MRR along the direction parallel to the best-fit line on 
the MR-plane and the direction perpendicular to it. The 
view angle is chosen in such a way that the projections on 
the plane parallel and perpendicular to the MRR line can 
be easily figured out. The projection on the MR-plane has 
already been displayed in Fig.[T] While the projection on the 
plane parallel to direction of the best-fit line can be easily 
understood in terms of selection effects (galaxy mass fall- 
off at the high mass end and lack of data at the low mass 
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Figure 2. 3D view of the MRR. The number frequency dis- 
tribution of galaxies on the MR-plane: M s is the total stellar 
mass in solar units and R1/2 i s the half-mass radius in kpc, 
and "Counts" is number of galaxies falling within each cell of 
the plane with dimensions AlogMi2=0.0 5 and A log _R 1 /? =(), 05 . 
The data are from the HB sample of [Bcrnardi ct al I l|2010l ). 
The thick line on the MR-plane is the linear fit of the data: 
logR 1/2 = 0.54 logAf s - 5.25. 

end), the distribution perpendicular to this is more difficult 
to explain. Indeed galaxies tend to fall in a rather narrow 
strip of the MR-plane, tightly gathering around the line with 
slope 0.54- 



3 NB-TSPH MODELS OF GALAXIES ON THE 
MR-PLANE 

Although a great deal of our analysis will be made using an- 
alytical relationships, some NB-TSPH models of ETGs by 
the authors are also considered. The models ar e presented in 
two g roups: the so-called reference models bv lMerlin et alj 
|2012f) which are calculated either to the present or down 
to redshift lower than z=l, and the ancillary models that 
are explicitly calculated for this paper and limited to ini- 
tial stage. In Appendix A we briefly report on a few salient 
properties characterizing all the models in use, whereas the 
complete description of the reference models and underly- 
ing e volutionary code can be found in iMerlin et all (2010, 
l2012h . 

For the sake of completeness, in this section we limit 
ourselves: 

(i) to remind the cosmological scenario adopted to 
calculate the galaxy models, i.e. the standard A-CDM 
model with H o =70.1 km/s/Mpc, flat geometry, n A =0.721, 
08=0.817 and the baryonic fraction ~ 0.1656. The initial 
conditions (positions and velocities of the DM and BM par- 
ticles) are taken from large scale cosmological simulations 
taking into account the expansion of Universe; 

(ii) to display in Fig.[3]the position of the reference and 
ancillary models on the MR-plane; 

(iii) to discuss two general aspects of the models that 
are relevant here, the filiation thread and the paths of the 
models in the MR-plane. 

The filiation thread. The models we have calculated at 
varying the initial cosmological density contrast 5pi(z) and 
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Figure 3. Theoretical MRR, i.e. log-Ri/2 versus logM s , where 
R\/2 is the half-mass radius and logM s the total stellar mass. 
The filled squares are the twelve reference models of Table [2] and 
the circles are the ancillary models of Tablel3lof Appendix A. The 
color code adopted for the various models is in principle related 
to the initial over-densities. Since the entries of Tables [2] and 
[3] are fully sufficient to identify the models, the correspondence 
color-code-initial density is superfluous here. The thick black line 
is the fit of the twelve reference models. The thin straight lines 
are the theoretical MRR on which proto-galaxies of given DM 
over-density settle down at the collapse stage as a func tion of the 
redshi ft. The lines are given by eqn. Q taken from iFan et al] 
(2010): solid lines are for m = 10 and / CT = 1; the color indicates 
different redshifts (black: z = 0; blue: 2 = 1; red: z = 2.5; green: 
z = 5; cyan: z = 20); the thin black dashed line is for m = 10, 
f a = 1.5, and z = 0; Finally, the thin black dashed-dotted line is 
for m = 5 and f a = 1, and z = 0. 



the star formation efficiency e s / help us to understand the 
different role played by 5pi(z) and the gas density p SlS / at 
the onset of star formation in determining the size of the 
resulting galaxy made of stars. Given the total mass oh the 
proto-halo Mdm, the cosmological density Spi(z) determines 
the initial radial dimension of the DM perturbation. This 
does not coincide with the initial radius at which star forma- 
tion begins in the baryonic component of a galaxy. In other 
words, within the potential well of DM, the gas keeps cooling 
at increasing density and only when the threshold density for 
star formation is met (which in turn is related to e s /), stars 
appear: the galaxy is detectable on the MR-plane. In this lat- 
ter step of the filiation thread, p 9 , s / plays the dominant role. 
The remaining gas continues to fall into the gravitational po- 
tential well until either it is exhausted by SF or it is expelled 
via shocks because of energy feed back. Thanks to this there 
will be a correlation between the DM and the final mass in 
stars, measured by m = Mdm/M s (wi th m ~ 10 on the 
average, see below). NB-TSPH models l|Merlin et al.ll2012T ) 
indicate that the transformation of BM into stars occurs un- 
der the homology condition GMdm/Rdm — GM S /R S , i.e. 
equal gravitational potential energy per unit mass of the two 
components. In general, the model galaxies follow this rule 
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and the two components of a galaxy lay on nearly parallel 
MRRs: i.e. Rom oc M D f M and R 3 oc Ml" 3 with /? ~ 3. At 
decreasing total mass, the exponent goes from 0.333 for 
galaxies with total mass 10 13 M Q to about 0.2 for a mass of 

or even lower. This deviation from the R s oc M s 
law can be interpreted as due to an increasing departure 
from the condition of an ideal collapse because dissipative 
processes are now at work. The higher the initial mass, the 
closer the evolution of the proto-galaxy is to the simple col- 
lapse models. In other words, the straight collapse configu- 
ration corresponds to a minimum total energy of the system, 
whereas in all other cases the total energy system is far from 
the minimum level. Recasting the concept in a different way, 
the straight collapse is favored with respect to other energy 
costing configurations. Real galaxies tend to follow the rule 
by as much as they can compatibly with their physical con- 
ditions (total mass, initial density, star formation history, 
...). The result of it will be that the final model galaxies will 
be located on a new line (not necessarily a straight line but 
likely close to it) with a certain slope (about 0.2 in our case, 
see the linear fit of the final mean location of the reference 
models). The slope is flatter than the 0.33 slope of the iso- 
density line of the initial DM proto-galaxies. There is not 
reason for the two slopes being the same0 

Path of the models in the MR-plane. Given these 
premises, a key question to clarify from the very beginning 
is the following: by how much our model galaxies, in which 
the so-called early hierarchical scheme is at work - early ag- 
gregation of lumps of DM and BM, conversion of gas into 
stars according to some star formation history, expulsion of 
gas by galactic wind over long periods of time - do change 
their star mass and dimensions during the whole evolution- 
ary history? To this aim we show in Fig. [6]the paths of some 
galaxy models per mass group (all the others have similar 
behavior). The circles of increasing size mark the stage at 
which the models have assembled 25%, 50% and 75% of the 
final stellar mass. Models move erratically in the MR-plane. 
The situation can be understood considering the definition 
of radius we have adopted: the galacto-centric distance in- 
side which half of the stellar mass in enclosed. The radius for 
the early evolutionary stages is rather uncertain. From an 
operational point of view, in principle at each evolutionary 
stage, once identified the barycenter of all the star particles 
existing at that time, spheres of increasing radius centered 
at this point are drawn up to encompassing half of the stel- 
lar mass. In the early evolutionary stages, characterized by 
captures of lumps of baryonic mass onto the already existing 
stellar body, desegregation of part of it by close encounters 
of other lumps of matter and similar events, the barycen- 
ter and the size of the sphere containing half star mass are 
quickly changing with time. The situation stabilizes as more 
and more baryonic mass is turned into stars and the main 
body of the galaxy is shaped. All this reflects on the the 
quite erratic path of the models on the MR-plane. However, 
looking at the large-scale trend, over the whole lifetime of 



2 In principle, the different redshift of the haloes at the collapse 
stage should be taken into account when considering the slope of 
their MRR. However, as our models have quite similar collapse 
redshifts, z ~ 2.5 — 5, the slopes of their MRRs do not change 
significantly. 
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Figure 4. Path on the MR-plane of three reference models of dif- 
ferent mass and same initail over-density. The open circles along 
the solid lines show the position at the time at which 25%, 50% 
and 75% of their stellar content is assembled; the size of the cir- 
cles is proportional to the assembled mass. As the stellar mass 
of the model galaxies grows monotonically with time, the path 
on the MR-plane is in general from left to right. However, some 
erratic displacements are possible owing to the definition itself of 
ill/2 > in particular during the earliest evolutionary stages. 

the order of 13 Gyr, the star mass can increase by a fac- 
tor of 10, whereas the radius remains nearly constant for 
the low and the intermediate mass galaxies. The situation is 
more confused for the high mass models. In any case, each 
model galaxy wonders around in a rather small box of the 
MR-plane. 

In this context we point out that, while e s f = 1 may 
concur to create diffuse stellar systems rather than compact 
objects because of the artificially high efficiency in periph- 
eral regions, the physical argument behind the result is ro- 
bust. The process of star formation most likely depends on 
the local properties of the medium, primarily in terms of gas 
density, which (at least as a first approximation) depends on 
the depth of the potential well. Therefore, a "threshold den- 
sity" (whether universal or not is still an open question) 
should be considered as an unavoidable ingredient in the 
scenario. Threshold density and efficiency parameter are un- 
known. However, they can vary without changing the phys- 
ical core of our conclusions: the position of a galaxy on the 
MR-plane is essentially decided in the early stages of its 
evolution. 



4 INTERPRETATION OF THE MR-PLANE 

If we compare the present-day position of the reference NB- 
TSPH models on the MR-plane with the region populated by 
real galaxies (Fig. [3j , at a first glance one would be tempted 
to conclude that only the high mass models fairly agree with 
observations, whereas the low mass ones (and to some extent 
also those of intermediate mass) apparently have too large 
radii with respect to their masses. However, before drawing 
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redshift 

Figur e 5. The growth function of haloes n(MoM , z) reproduced 
from lLukic et alJ l|2007l l. 



tion is expressed in terms of the normalized Hubble con- 
stant h = //o/100, where Ho is assumed to be -Ho=70.1 
Km/s/Mpc. The explored inter val of redshift goes from to 
20. The n(M DM , z) function of iLukic et all (|2007l ) is shown 
in Fig. Of. 

Although w hat we are going to say is w ell kn own, see the 
pionee r study of IPress fc SchechteJ [|l974h and ILukic et alJ 
(2007, for ample referencing), for the sake of clarity and as 
relevant to our discussion we note the following: (i) for each 
halo mass (or mass interval) the number density is small at 
high redshift, increases to high values toward the present, 
and depending on the halo mass either gets a maximum 
value at a certain redshift followed by a decrease (typical 
of low mass haloes) or it keeps increasing as in the case of 
high mass haloes; in other words, first creation of haloes of 
a given mass (by spontaneous growth of perturbation to the 
collapse regime or by mergers) overwhelms their destruction 
(by mergers), whereas the opposite occurs past a certain 
value of the redshift, for low mass halos; (ii) at any redshift 
high mass haloes are orders of magnitude less frequent than 
the low mass ones; (iii) at any redshift, the mass distribution 
of haloes has a typical interval of existence whose upper mass 
end (cut-off mass) increases at decreasing redshift. 



the conclusion that essentially the models fail to reproduce 
the data, it is worth recalling that similar trend is first pre- 
dicted by the filiation thread, second it has bee n found by 
other NB-T SPH calculations, IChiosi fc Carrarol (see for in- 
stance Hqq3). 

Is the observational sequence populated only by galax- 
ies behaving as our massive ones? Or what else? Since for 
any value of the halo mass there is a certain redshift be- 
low which haloes of this mas s start decreasing in number 
by mergers ijLukic et al.ll2007h . galaxies generated by those 
haloes become more and more unlikely as it should be the 
case for our low mass models. Indeed when at a given red- 
shift we have assumed the existence of haloes of any mass, 
we have neglected this important effect. Similar considera- 
tions would also apply to haloes of large mass. Therefore, 
the situation may occur that haloes/baryonic galaxies are 
calculated and plotted onto the MR-plane even though ac- 
cording to the above arguments their existence is very un- 
likely. On this ground, in the following we argue that the 
observational MRR of ETGs (galaxies in general) is the re- 
sult of convolving two agents: the halo growth function pro- 
viding the number density of haloes of different mass as a 
function of the redshift (the concordance A-CDM Universe), 
and the fundamental MRR determining the size of a galaxy 
as a function of its mass and formation redshift. These two 
important ingredients are shortly presented below. 



4.2 The fundamental Mass-Radius Relationship 

Spherical DM perturbations (with BM inside) that undergo 
collapse when the density contrast with respect to the sur- 
rounding medium reaches a suitable value obey the MRR 
given by 

where p u (z) is the density of the Universe at the redshift z, 
and A is the factor for the density contrast of the DM halo. 
This expression is of general validity whereas the function A 
depends on the cosmological model of the Universe, includ- 
ing the A-C DM case. All details and demonstration of it can 

be found in (|Brvan fc Normanlll99ol their Eq. 6). 

In the context of the A-CDM cosmology, iFan et all 

(2010) have adapted the general relation (J5J to provide an 
expression correlating the halo mass Mo m and the star mass 
M s of the galaxy born inside it, the half light (mass) radius 
Ri/2 of the stellar component, the redshift at which the col- 
lapse takes place Zf, the shape of the BM galaxy via a co- 
efficient Ss{ns) related to the Sersic brightness profile from 
which the half-light radius is inferred and the Sersic index 
ns, the velocity dispersion of the BM component with re- 
spect to the that of DM (expressed by the parameter / CT ), 
and finally the ratio m = Mom /Ms- The expression is 



4.1 The halo growth function u(Mom, z) 

The distribution of the DM halo masses and their rel- 
ative number density as a function of the redshift has 
been recently studied, among the others, by ILukic et al.l 
ll2007h who, using the A-CDM cosmological scenario and 
the lWarren et all(|2006t ) mass function of haloes, derive the 
halo growth function ti(Mdm , z). This gives the number 
density of haloes of different masses per (Mpc//i) 3 result- 
ing by all creation/destruction events. The growth func- 



o _ nQ Ss(n)25 /1-5\V M DM V /3 4 

Typical value for the coefficient Ss(ns) is 0.34. For the ratio 
m = Mom /M s the empirical data confine it in the range 20 

3 This has been derived fr om an analy t ical in terpolation of the 
data presented in Fig. 1 of ILukic et al ] l|2007l l. More details arc 
given below. 
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to 40, whereas our NB-TSP models yield m ~ 15 to 20 de- 
pending on the particular case we refer to. For the purposes 
of the present study we round the ratio to 10 (see below). Fi- 
nally fcr yields the three dimensional star velocity dispersion 
as a function of the DM velocity dispersi on, a s = f a apM - 
Here we adopt f a = 1. For more details see iFan et all (|2010h 
and references therein. 

We point out that this relation is strictly valid only 
for monolithic infall of BM into collapsing DM potential 
wells. Therefore, its application is not justified when dealing 
with late mergers between already virialized objects, or with 
different mechanisms of mass assembly, such as, in particu- 
lar, an early-hierarchical scenario. Nevertheless, this formula 
provides a general reference to obtain the typical dimension 
of a galactic system as a function of its mass and formation 
redshift. While adjustments are possible, the general trend 
is well defined. However, some deviations from this law are 
possible and expected, e.g. for low redshifts. See below for 
further discussion. 

The slope of relation Q is nearly identical to (only 
slightly steeper than) the slope estimated from the NB- 
TSPH models; the difference can be fully ascribed to the 
complex baryon physics, which causes the stellar system to 
be slightly offset with respect to the locus analytically pre- 
dicted from DM haloes. Therefore, a model slope (close to 
1/3) different from that of the observational MRR is not 
the result of inaccurate description of the physical processes 
taking place in a galaxy; on the contrary, it mirrors the 
fundamental relationship between mass and radius in any 
system of given mean density. Indeed it is remarkable that 
quite complicated numerical calculations clearly display this 
fundamental feature. If this is the case, why do real galaxies 
gather along a line with a different slope? 

4.3 The cosmic galaxy shepherd 

To answer the above question, we attack the problem from a 
different perspective, trying to investigate whether the ob- 
servational MRR owes its origin to deeper re asons, likely 
relate d to the growth function of DM haloes l)Lukic et al.l 
l2007f) . 

Given a certain fixed number density of haloes 
N s (thereinafter referred to as a "statistics"), on the 
n{Mr)M, z) — z plane of Fig. [5] this would correspond to an 
horizontal line intersecting the curves for the various masses 
at different redshifts, i.e. obeying the equation n(MoM, z) = 
N 3 . Each intersection provides a pair (Mdm, z) which gives 
the mass of the haloes fulfilling the statistics N 3 at the cor- 
responding redshift z (or viceversa the redshift satisfying 
the statistics for each halo mass). For any value N s we get 
an array of pairs (Mdm, z) that can b e extrapolated to 
a continuous function that, with aid the IFan et all (|2010l ) 
relationship (in which the parameters m and f a are fixed), 
provides the corresponding relationship between the mass in 
stars and the half-mass radius of the baryonic galaxy asso- 
ciated to a generic host halo to be plotted on the MR-plane. 

Repeating the procedure for different values of N s , we 
get a manifold of curves on the MR-plane. It turns out 
that with the "statistics" corresponding to 10~ 2 haloes per 
(Mpc//i) 3 , the curve is just at the edge of the observed distri- 
bution. Higher values of the statistics would shift it to larger 
haloes (baryonic galaxies), the opposite for lower values of 




8 9 10 11 12 13 
log(M s ) 

Figure 6. The cosmic galaxy shepherd (CGS) and the corre- 
sponding locus of DM parent haloes (the red thick and thin solid 
lines, respectively). The "statistics" is 10 — 2 haloes per (Mpc/h) 3 . 
In addition to this we show the case with "statistics" of 10~ 8 
haloes per (Mpc/h) 3 (the magenta thick and thin dashed lines). 
The various loci are plotted onto the observational MR-plane (see 
text for all det ails). We also d r aw th e observational data for the 
HB sample of iBernardi et al] ll2010h with their linea r fit (solid 
black line), and two theoretical MRR from eqn. {4} bv lFan et alj 
(2010) (solid thin blue lines), relative to z = and z = 10, with 
m = 10 and fa- = 1. The present day position of the reference 
galaxy models and their linear fit are shown (filled squares and 
black solid line). Finally, the filled black circles show the inter- 
section of the CGS with the theoretical MRR for the case with 
10~ 2 haloes per (Mpc//i) 3 . 

the "statistics". Why is N a = 10~ 2 haloes per (Mpc//i) 3 so 
special? Basing on crude, simple-minded arguments we re- 
call that the total number of galaxies observed by the SDSS 
amounts to about ~ 10 6 , whereas the volume of Universe 
covered by it is about ~ 1/4 of the whole sky times a depth 
of ~ 1.5 x 10 9 light years, i.e. ~ 10 8 Mpc 3 , to which the 
number density of about 10~ 2 haloes per (Mpc//i) 3 would 
correspond^. 

The analytical fit of this curve is 

log.Rr/2 =0.048562(logM s ) 3 - 1.4329(log M s f 

+ 14.544(log M s ) - 50.898 (5) 

and we note that the slope gradually changes from 0.5 to 
1 and above as we move from the low mass to the high 
mass range. It is worth recalling here that a similar trend 
for the slope is also indic ated by the observational data (see 
Ivan Dokkum et al.ll201ol . and references therein). Owing to 

4 We are well aware that this is a very crude estimate not tak- 
ing into account many selection effects both in the observations 
and the ha l o stat istics based on NB simulations, such as the 
iLukic et al.l (I2007T) plane itself. However, just for the sake of ar- 
gument, we can consider it as a good estimate to start with. 



8 C. Chiosi, E. Merlin and L. Piovan 



the many uncertainties we do not try to formally fit the me- 
dian of the empirical MRR, but we limit ourselves to show 
that the locus predicted by the 10~ 2 haloes per (Mpc/fe) 3 
"statistics" falls on the MR-plane close to the observational 
MRR. Lower or higher values of the halo number density 
would predict loci in the MR-plane too far from the obser- 
vational MRR unless other parameters of eqn. Q, e.g. m, 
are drastically changed assuming values that are difficult to 
justify. 

Finally, we call attention on the fact the locus on the 
MR-plane defined by the relation ([5j is ultimately related 
to the top end of the mass scale of haloes (and their fil- 
iated baryonic objects) that can exist at each redshift. In 
other words, recalling that the mass of any intersection pair 
for N s — 10 -2 corresponds to haloes becoming statistically 
significant in number on the observed spatial scale at the as- 
sociated redshif t, this can be interp r eted a s the so-called cut- 
off mass in t helPress fc Schechtej (| 19741 ) or equivalent for- 
malisms fsee lLukic et al. 20071 . for details and references). 
Therefore, this provides also an upper boundary to the mass 
of galaxies that are allowed to be in place (to collapse) at 
each redshift. We name these locus the Cosmic Galaxy Shep- 
herd (hereafter CGS). All this is shown in Fig. [6] where we 
also plot the curves relative to another possible "statistics" 
- that is, 10 -8 haloes per (Mpc//i) 3 , corresponding to 1 halo 
per 10 s (Mpc//i) 3 , for the sake of comparison. 

There are two points to be clarified. First, this way of 
proceeding implies that each halo hosts one and only one 
galaxy and that this galaxy i s an early type o bject matching 
the selection criteria of the \Bernardi et all \20lA) sample. 
In reality ETGs are often seen in clusters and/or groups of 
galaxies and many large spirals are present. Only a fraction 
of the total population are ETGs. One could try to correct 
for this issue by introducing some empirical statistics about 
the percentage of ETGs among all types of galaxy. Despite 
these considerations, to keep the problem simple we ignore 
all this and stand on the minimal assumption that each DM 
halo host at least one baryonic component made of stars. 
This is a strong assumption, on which we will come back 
again later on. Second, we have assumed m = 10 and / CT 
I , According to iFan et al.l i|2010l '). the empirical estimate of 
Mdm /M a ration is about 20-40, i.e. a factor of two to four 
less efficient star formation than we have assumed basing on 
our NB-TSPH models. However, a smaller value for m does 
not invalidate our analysis, because it would simply shift 
the location of the baryonic component on the MR-plane 
corresponding to a given halo "statistics". Finally, / CT = 1 
is a conservative choice. The same considerations made for 
m apply also to this parameter. At present, there is no need 
for other values. 

Along the CGS, redshift and cut-off mass go in inverse 
order, i.e. low masses (and hence small radii) at high 
redshift and viceversa. This means that a manifold of 
MRRs defined by eqn. @, each of which referring to a 
different collapse redshift, can be selected, and along each 
MRR only masses (both parent Mdm and daughter M 3 ) 
smaller than the top end are permitted, however each of 
which with a different occurrence probability: low mass 
haloes are always more common than the high mass ones. 
In the observational data, it looks as if ETGs should occur 
only towards the high mass end of each MRR, i.e. along 
the locus on the MR-plane whose right hand side is limited 



by the CGS. This could be the result of selection effects, 
i.e. (i) galaxies appear as ETGs only in a certain interval of 
mass and dimension and outside this interval they appear 
as objects of different type (spirals, irregulars, dwarfs etc.), 
or (ii) they cannot even form or be detected (e.g. very ex- 
tended objects of moderate/low mass). Finally, in addition 
to this, we argue that another physical reason limits the 
domain of galaxy occurrence also on the side of the low 
mass, small dimension objects. We will come to this later on. 

Dissipation-less Collapse. It is an easy matter to under- 
stand that the CGS is another way of rephrasing the top-hat 
sp herical dissipa t ion-le ss collapse for primordial fluctuations 
bv lGott fc Reel jHzl). For the sake of easy understanding 
we summarise the key ste ps of t h e diss ipati on-less collapse 
using the version given bv lFaberi dl98i l and lBurstein et al.l 
(1997). Let 8 be the rms amplitude of primordial density 
perturbations 



8 oc M"5" 



(6) 



where M is the mass at the initial red-shift, and n is the slope 
of the density fluctuation 8. After collapse, the equilibrium 
structure of a DM halo originated from given S and M follows 
the relations |Gott fc Reesll 19751 ) 



RtxS^M* 
from which we immediately get 



R oc m"^" 



(7) 



(8) 



As already pointed out bv lBurstein et all (Il997), inserting 
n = —1.8, the power spectrum of CDM (|Blumenthal et all 
Il984h . we get the relation 



Rdm oc M° dm 



(9) 



The slope of the MRR derived from the dissipation-less col- 
lapse is the same of eqn. ([TJ all over the mass range from nor- 
mal/giant galaxies, M32 and o;Cen like objects, to classical 
Globular Clusters . This was already pointe d out long ago b y 
Chiosi fc Carrarol (|2002T ) and discussed bv lGrahaml i|201ll ). 
It can be easily checked by plotting the data for Globular 
Clusters, M32 and ojCen in Fig. [T] The advantage of the 
CGS with respect to the simple dissipation-less collapse is 
that it provides slope and zero point of the observational 
MRR and also predicts its change in slope at increasing star 
mass of the galaxy. 

Given these premises, we suggest that the observational 
MRR, eqn. JTJ, represents the locus on the MR-plane of 
galaxies whose formation and evolution closely followed the 
scheme of dissipation-less collapse, i.e. the ones of the largest 
mass for each formation redshift. The dwarf galaxies have 
a different interpretation, because they significantly depart 
from the above evolutionary scheme and the MRR holding 
for ETGs. As already pointed out they follow a much flatter 
relationship which mimics the theoretical MRR (slope 1/3) 
even if this may be a mere coincidence. 
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Table 1. Coefficients of the polynomial interpolation of the relation l|10jl , which provides the number density of haloes n(MoMi z ) P er 
(Mpc/h) 3 . 



Mass [M@/h] 


A 4 


A 3 


A 2 


Ai 


Ao 


5e7 


-2.34275e-5 


1.28686e-3 


-2.97961c-2 


2.11295e-l 


2.02908 


5e8 


-2.76999e-5 


1.49291e-3 


-3.47013c-2 


2.13274e-l 


1.13553 


5e9 


-1.31118e-5 


6.50876e-4 


-2.36972e-2 


1.31993e-l 


0.23807 


5el0 


-1.18729e-5 


6.65488e-4 


-3.17079e-2 


1.30360e-l 


-0.59744 


5ell 


-1.47246e-5 


8.10097e-4 


-4.65279e-2 


1.13790e-l 


-1.44571 


5el2 


6.59657e-5 


-7.19134e-4 


-6.99445e-2 


1.06782e-l 


-2.45684 


5el3 


-7.34568e-4 


9.99022e-3 


-1.65888e-l 


-9.48292e-2 


-3.11701 


5el4 


4.89975e-3 


-5.17004e-2 


-1.61508c-l 


-5.83065c-l 


-4.28270 



5 SIMULATIONS OF THE MR-PLANE 

The CGS suggests that a deep relation exists between the 
way galaxy populate the MR-plane and the cosmological 
growth of DM haloes and furthermore it accounts for the 
slope of the observational MRR. However, it does not pro- 
vide an explanation for the tightness of the MRR. To cast 
light on this issue, we analyze the whole MR-plane, apply- 
ing the same line of reasoning used above, and present a 
simulation of the MR-plane and the "observational" MRR 
derived from first pr i nciple s. Once again, we start from the 
study of iLukic et all (|2007h and fit their n(MoM , z) curves 
with fourth order polynomial expression^]. Then we count 
the total number of haloes per mass-bin A log Mo m at red- 
shift z = 0. This is simply given by reading off the val- 
ues of the curves along the y-axis and interpolating for in- 
termediate values (we take equally spaced logarithmic bins 
AlogAf = 0.05). These are the haloes that would nowadays 
populate the synthetic MR-plane and that should be com- 
pared with the observed galaxies. Of course, to make a mean- 
ingful comparison , as the total number of haloes read off the 
ILukic et all (|2007h plot refers to a volume of 1 (Mpc//i) , one 
has to scale it by a suitable factor C to match the real vol- 
ume of the observed portion of the sky from which the data 
are obtained. We estimate that C = 5 x 10 6 is a reasonable 
choice. This is based on the following arguments: the ob- 
servational sample we use contains about ~ 60, 000 objects, 
i.e. 6% of the total number of objects observed by SDSS 
(~ one million); therefore, we should consider ~ 5 — 10% 
of the total number of haloes expected in the total volume 
of the sky observed in the survey, which turns out to be 
~ 10 s Mpc 3 , so C ~ 0.05 x 10 s . It is worth recalling here 
that procedure is the same as the one we have been using 
to determine the statistics N B = 10 -2 generating the CGS 
corresponding the observational MRR. Once more, the line 
of reasoning we follow may seem very crude, but it is fully 
adequate to the purposes. 

Now we have to assign each halo of mass Mom a radius 
Ri/2 using eqn. @. The radius is also function of the for- 
mation redshift Zf. Therefore, we must assign a typical for- 
mation redshift to each halo of mass Mom to be able to plot 



5 I ndeed this p o lynom ial fit has already been used to reproduce 
the iLukifTe t al. (2007j) plane shown in Fig. [5] 



the corresponding visible galaxy on the MR-plane. To this 
aim, we pro ceed as follows. Firs t, we start from the growth 
functions of ILukic et all |200?tl and derive the number of 
tl{Mom,z) of haloes existing in each mass bin Alog Mom 
at any redshift z. This number is the result of two compet- 
ing effects: the formation of new haloes of mass Mom via 
merger and/or acquisition of lower mass haloes, and the de- 
struction of haloes of mass Mo m because of they merge to 
form higher mass haloes. Therefore, the following equation 
can be written 



n(M D M, z) =n(M D M, z + Az)+ 

n+(M D M,z) -n-(M DM ,z + Az) (11) 

where n+ and n_ represent the creation/destruction mech- 
anisms. In particular, the quantity we are interested in is 
n+(MoM , z), which is the number of new haloes of mass 
Mom which are born at redshift z. 

The number of haloes that merge to form higher mass 
systems is in turn a fraction of the number of haloes existing 
at that time, i.e. U-(Mom , z + Az) = rj x u(Mom, z + Az), 
with < rj < 1; so 

n(M D M, z) =n(M D M, z + Az) + 

n + {M DM , z)-rjX n(M D M, z + Az) (12) 

and 

n + {M D M,z) =n(M D M, z)- 

(1 - 77) x n(M D M, z + Az) (13) 

The only free parameter here is rj, the fraction of haloes 
that merge to form higher mass systems in the redshift in- 
terval Az — 0.145 x z + 0.1 0. In principle, the fraction rj 
could vary with the redshift. However, for the sake of sim- 
plicity we assume that rj remains constant. Thus, we obtain 
a value N+(Mdm , z) for each interval Mom , Mom + AMom 
and z, z + Az. This number, re-normalized to unit over the 
whole interval, can be considered as the relative probability 
that a halo of mass Mom is born at redshift z. 

Finally, for each halo of mass Mom we compare a ran- 
domly chosen number q G [0, 1] with the cumulative proba- 
bility 



n(M DM ,z) = A n (MoM) x z'\ 
The coefficients A„(Mom) are listed in Table l4~3l 



(10) 

6 The mass and redshift intervals are as in ILukic et al. I l l2007l) . 
The expression for Az secures that nearly equally spaced intervals 
are used for logi?. 
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z ~ z i 

P ^ = n + (M DM ,z) 

z — z max 

until we have q < P Zi , and take zf — Zi as its formation red- 
shiffl According to eqn. Q, once m and f a are known (we 
adopt here 10 and 1, respectively), the radius of the bary- 
onic galaxy filiated by this halo, R\/2{Mdm, zf), is known 
and the galaxy can be plotted on the MR-plane. 

Given these premises, we derive the galaxy counts to 
be compared with with the observational data. To this aim, 
we divide the MR-plane in a large number of elemental cells 
with dimensions AlogM s = Alogi?!/2 = 0.05, assume a 
value for the free parameter 77, calculate the number of galax- 
ies expected to fall in each cell, and finally compare this with 
the observational counterpart (our SDSS sample). 

Before looking at the results, there is an important issue 
to clarify: the distribution of galaxies on the MR-plane is 
found to heavily depend on the fraction 77 of haloes which 
merge to form higher mass systems in the redshift interval 
Az. The point has been checked by adopting different values 
of r\ and comparing the results with the observational data. 
We find that in general the lower the value ofrj, the better the 
simulations agree with the observational data. For instance, 
with the present spacing of the redshift, 77 = 0.01 yields 
excellent results. The agreement would be even better for 
77 = 10 -4 . In contrast for larger values of 77, e.g. 77 = 0.05 (5% 
of the haloes in each mass bin merge during each redshift 
interval) the results significantly worsen. In the following, we 
will limit to the case with r\ = 0.01 (1% of the haloes merge 
during each redshift interval Az), whereas the case with 77 = 
0.05 and an attempt to explain this trend is postponed to 
Section T5. 21 below. The case 77 = 0.01 can be considered as a 
sort of acceptable upper limit for this parameter. 

The results for the case 77 = 0.01 are plotted in the 3D 
space log M s , log R1/2, n(M 3 , R1/2) (otherwise referred to as 
"Counts") of Fig. [7] in which the previous variables Mdm 
and z are now replaced by M„ and Ri/2- Along the verti- 
cal axis, the positive values (blue points) show the observa- 
tional data, whereas the negative values show the theoretical 
counts (red points). In this way, the one-to-one comparison 
is straightforward. The view angle is chosen in such a way 
that the peak values, fall-off toward the high mass range, 
and thickness of the MRR are simultaneously visible. 

We do not make use of the low mass end of the MRR 
because the observational one is obviously affected by incom- 
pleteness and other selection effects. Indeed, while the the- 
oretical expectation is that the number of low mass haloes 
(low mass visible galaxies) should increase toward the low 
mass end, the opposite is seen in the data. This difference is 
most likely due to observational incompleteness even though 
other physical causes inhibiting either the formation of stars 
in low mass haloes (thus making the galaxy hardly observ- 



7 The upper limit for the redshift is not always z max = 20, but 
the value given by the relation 

z max = INT(0.0126M 3 - 0.7597M 2 + 3.5848M + 17.571) (14) 

with M z = Mi — (M m i n — 1) obtained by fitting the lLukic et alj 
1120071 1 data. 
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Figure 7. Comparison of the theoretical (red dots, negative z- 
axis) with the observational (blue dots, positive z-axis) number 
frequency distributions of galaxies in the MR-plane: log M5 is 
the total stellar mass (in solar units) and log(i?!/2) the half- 
mass radius (in kpc), "Counts" is the quantity An(z) = n(z + 
Az) — n(z), i.e. the number of haloes per Mpc 3 that are born 
between 2 and 2 + Az plotted per each elemental cell of the MR- 
plane. The thick line is the best fit of the observational data: 
logR^j2 = 0.54 l ogMg — 5.25. The ob servational data are from 
the HB sample of Bernardi ct al ] lj20ld l. The theoretical number 
frequencies have been plotted as negative quantities to get the 
mirror image of the observational data. 

able) or even destroying part of the low mass haloes cannot 
be excluded. Related to this, the discrepancy between pre- 
dicted and observed numbers of dwarf galaxies has long been 
debated with no firm conclusions. Current hypotheses sug- 
gest that small haloes do exist, but fail in producing stars 
and luminous matter. See below for more details on this 
issue. 

Furthermore, we note that the high-mass tail of the 
theoretical MRR is too populated with respect to the ob- 
servational data. This is likely to be a consequence of the 
assumption that one halo always hosts one and only one 
galaxy. In the real Universe, massive haloes may actually 
host many galaxies, because merging haloes does not im- 
ply merging of their galaxy content. For example, galaxy 
clusters (of typical masses M c i U3ter — 10 14 — lO 15 A/0) are 
known to host hundreds or thousands of galaxies; groups of 
galaxies can host tens of galaxies, some of which may share 
a unique dark halo. Thus, in the high-mass tail of our sim- 
ulation large galaxies corresponding to massive haloes may 
be replaced by smaller, less massive galaxies, which would 
reconcile theory with data. No attempt is taken to cure the 
theoretical MRR. In any case, the similarity between real 
and simulated counts is striking. 

Looking now at the the thickness of the theoretical and 
observational MRRs, at any given radius, while on the high 
mass side the MRR is bounded by the CCS, its low mass side 
is determined by the fact that at each halo mass there is a 
certain redshift below which the number density n(MoM, z) 
starts decreasing: low mass haloes decrease because they ei- 
ther merge into bigger objets or are captured by more mas- 
sive haloes. As a consequence, the MRR is expected to be 
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a strip on the MR-plane, populated by galaxies of different 
type, among which toward the high mass end ETGs (whose 
structure and evolution closely resembles the dissipation-less 
collapse) are more numerous. This is clearly shown by the 
simulation. 

After this analysis of the whole log[M 3 ], logfT?^], 
N(Ms, R1/2) space, we simulate the MRR of galaxies by 
drawing as many points in each cell of the MR-plane as in- 
dicated by the Counts displayed in Fig. [71 we blur the results 
by artificially applying a small random displacement to each 
point, considering possible small variations in the m and / CT 
parameters. The results are plotted in Fig. [5] (left panel), 
and c ompared to the real MR-plane of the iBernardi et al.l 
(2010) sample (right panel). Again, the agreement is good; 
first of all, the simulated points fall exactly in the same re- 
gion in which real galaxies are found; second, the shape of 
the two distributions is similar; third, the two distributions 
have similar width. 

As expected, the theoretical distribution is broader than 
the observational one. The cause is merely of statistical na- 
ture. Indeed there is still something missing so that the sim- 
ulated MRR might be the real counterpart of the obser- 
vational data. In brief, the previous analysis has estimated 
the probability that a cell of the MR-plane corresponds to 
a galaxy with that mass and that radius. To get a realistic 
simulation of the MRR one has to translate the occupation 
probability of a cell into a number of real galaxies in the 
same cell. This is particularly true for those cells in which 
the probability is very low. While on the theoretical side 
these cells may have a small yet finite occupation probabil- 
ity, in reality with finite numbers of galaxies these cells may 
be empty. So in reality the MRR is expected to be thinner 
than predicted by theory. In order to take this into account, 
we need Monte-Carlo simulations of the MRR. This will be 
the subject of Section [5.11 below. 

Finally, the simulated MRR extends towards the realm 
of the low mass objects and does not find an observational 
counterpart in the data to dispo sal. This is less of a pr oblem 
as the observational sample by IBernardi et al.l (j2010) con- 
tains only objects brighter than a threshold magnitude, so 
that low mass (and hence low luminosity) galaxies are simply 
missing. Indeed we have data for dwarf galaxies only in the 
Local Universe, because observing dwarf galaxies in remote 
regions of the Universe is still out of reach. Thus, we may 
argue that many low mass objects should exist, but they are 
not detected, and hence they are not present in the observed 
sample. Moreover, we must also keep in mind that our simu- 
lations of the MRR stand on the assumption that each DM 
halo hosts a luminous baryonic object. While for intermedi- 
ate and high (but not too high!) masses this sounds reason- 
able, at the scale of low mass small galaxies the assumption 
may brake down. Many low mass haloes could indeed not be 
able to form visible galaxies, for instance if the primordial 
gas cannot collapse into the potential wells of these haloes, 
because after re-ionization the gas pressure remained high. 
In other cases, only diffuse, low surface brightness stellar 
systems may form. Finally, a number of physical processes 
may concur to dislocate a low mass galaxy from region of 
the MR-plane in which it was born to other regions of it: 
for instance, galactic winds, gas stripping, and others (see 
below) . 

However, we point out again that the number of small 
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Figure 8. Loft panel: Predicted MRR limited to galaxies more 
massive than 10 10 Mq fo r the case with t? = .01. Right panel: 
The observational data of IBernardi et all feoioh 



objects which end up populating the simulated MR-plane 
is consistent with the expectations from a hierarchical theo- 
retical cosmology; indeed, considering the ratios in the num- 
ber of haloes of diff erent masses at z = predicted by the 
Warren et all l|2006l ) mass function, the agreement with our 
simulation is excellent. Moreover, the same ratios are consis- 
tent with those between the number of luminous objects in 
the Local Group, considering the two massive galaxies, the 
dwarfs, and estimating the number of Globular Clusters. 



5.1 Monte-Carlo make-up of the MRR 

To get a simulation as close to reality as possible, we make 
use of the Monte-Carlo technique and consider u+(Mdm, z) 
as the probability that an halo (and the associated bary- 
onic galaxy) with the chosen mass and corresponding radius 
actually comes to life. 

To this aim, we pick a set of S random numbers r £ 
[0, 1], and compare each of them with the value P — C x 
ti+(Mdm, z) relative to a randomly chosen cell in the MR 
plane (C being a suitable normalization factor; the cells are 
defined as described in Sec. [5}. If r < P, a point is plotted 
on the plane at the position of the relative cell (with a small 
random displacement to blur the result), otherwise it is not. 
The result, adopting S = 3 x 10 5 and C = 10" 4 , is plotted 
in Fig. [5] together with the usual observational sample and 
with a few theoretical relations. 

According to our simulations of the MRR, there should 
be galaxies born at different redshift and hence objects of dif- 
ferent age. Fig. [TO] we group the galaxies according to their 
formation redshift. The mean age in each strip increases 
moving from the upper to the lower bou ndary of the MRR. 
This a grees with the trend found by IValentinuzzi et all 
(|2010al fbh in galaxies of the WING survey of X-selected 
galaxy clusters. 
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Figure 9. Monte-Carlo simulation of the MRR: comparison of 
the theoretical (black dots) with the observational (red dots) dis- 
tributions of galaxies in the MR-plane log(i?j/2) m kpc versus 
logAf s in units of Mq. The theoretical simulation is for rj = 0.01 
(see the text for details). Superposed to it are the reference models 
(filled squares) for intermediate and high mass galaxies and their 
linear fit, the CGS for the 10 — 2 objects per (Mpc/h) 3 "statistics": 
the thick dashed green line for the stellar component and the blue 
lines for the Da rk Matter halo. The thin trasverse lines are the 
I Fan et al. relationships for different redshift, namely z=0, 

5, 10, and 20 from the top to th e bottom. Fina l ly, the observa- 
tional data are the HB sample of lBernardi et al] fcoidl ). 




log[M s ] 

Figure 10. The theoretical MRR, log(R 1 / 2 ) in kpc versus log M s 
in Mq; the color code corresponds to different formation redshifts: 
cyan z < 0.5, green 0.5 < Z < 1, magenta 1 < z < 2.5, blue 
2.5 < z < 5, and black z > 5. The mean age of galaxies increases 
from the top to the bottom. 
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Figure 11. The same as in Fig. [9] but for rj = 0.05. 



5.2 Merger percentage 

The agreement we have reached so far has been possible 
using rj = 0.01. As already anticipated, while there is no 
problem with smaller values of rj, severe difficulties arise in- 
creasing it above the value we have adopted. For the sake of 
argument we show the results for rj = 0.05. The least we can 
say is that theory and data simply disagree. The qualita- 
tive explanation of this counter-intuitive result is as follows. 
Increasing rj (while keeping constant the total number of 
haloes of given mass Mom and redshift z) induces an in- 
crease of the number of haloes of the same mass Mom that 
collapse exactly at the same redshift z; in other words, if 57 is 
high many haloes merge to form even more massive haloes. 
Therefore, most of the haloes of mass Mom must be newly 
collapsed structures. On the contrary, low values of rj imply 
that fewer haloes disappear because of mergers, so that a 
larger fraction of the haloes in place at redshift z must have 
formed at earlier times and must have survived up to the 
epoch z. As a consequence of this complex game, the distri- 
bution function at the epoch of formation of halos with mass 
Mom changes shape, favouring more recent formation times 
as rj increases. In turn, recent formation times imply larger 
radii (eqn. 3}, thus displacing the visible galaxies toward 
the upper regions of the MR-plane. Basing on these results 
(excellent agreement between theory and observational data 
only for rather small values of rj), we are tempted to con- 
clude that mergers cannot have plaid a crucial role in the 
galaxy formation mechanism. 

To cast light on this issue, we assume that rj is constant 
with time and calculate the ratio between the total number 
of merger events, N merg — J2 Z . rj x ti(Mdm, z), and the 
total number of galaxies that ever existed in the Universe 
over the Hubble time, Ntot = Ylz n+(MoM , z), 
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(15) 



Crude calculations using the lLukic et al.l |2007h curves yield 



N„ 



1 /Ntot 



0.10 - 0.15 for rj = 0.01. This percentage 



increases to about 50% for rj = 0.05. With the latter value, 
the price to pay is that the corresponding MRR is much 
more dispersed than the observational one. Therefore, our 
analysis seems to suggest that only a modest fraction of 
haloes should have merged during the whole history of the 
Universe. Owing to the many uncertainties still affecting the 
above discussion, we do not insist on this issue. However, it 
is an important result of this study deserving a thorough 

investigation. 

Different conclusions are reached by iFakhouri et alj 

(|201Ct) who, using the joint data from the Millennium 
and Millennium-II simulations, reconstruct the past merger 
histories of haloes of different mass and estimate their 
merger frequencies. For halo masses in the rang e 10 10 M to 
10 15 Mq, ratio of progenitor mass in the range 10 -5 < £ < 1, 
and redshift interval < z < 15, they find the merger rate 
per unit redshift to be nearly independent of z and halo 
mass, and to depend on £ as oc £~ 2 . Finally, referring to 
a typical halo mass of 10 12 M ot j ot today, they estimate the 
probability of a major merger since z—1, 2 and 3 to be about 
31%, 53% and 69%, respectively. Haloes with the same mass 
at z—2, have about 50% probability of major mergers since 
z=4. The above probabilities refer to a specific halo mass, 
so the translation to the percentage of mergers with respect 
to the whole galaxy content over the entire history of the 
Universe is not straightforward. 

The only remark we can make, is that perhaps the two 
estimates could indeed be closer than it might seem at a 
first sight. The merit of our estimate is that observational 
data are used to compare the model r esults (choice of neon- 
strained on the MRR) , whereas the IFakhouri et al.l (|2010| ) 
analysis and estimates of the merger percentage are theo- 
retical, with little reference to the real world of galaxies we 
see in the Universe. Surely the topic is very alive and worth 
being pursued. 



5.3 Possible effects blurring the MRR 

We have shown that a tight MRR may result from first 
physical principles coupled to elementary cosmology. For the 
sake of completeness, we may turn the argument around and 
check whether starting from an intrinsically sparse distribu- 
tion of galaxies on the MR-plane, one may recover a tight 
MRR as the result of some physical effects. In other words, 
there are at least two important phenomena at work, i.e. 
mergers of two galaxies of different mass (and size) into a 
single object and matter ejection by a single galaxy by sud- 
den injection of energy (e.g. AGN activity and/or galactic 
winds) that could change the position of a given galaxy on 
the MR-plane and therefore scatter the expected distribu- 
tion. Whether, on the contrary, starting from a scattered 
distribution of galaxies, they can give rise to more or less 
ordered sequence such as that indicated by the MRR is hard 
to prove, and unlikely to occur in the real world. 

Minor and major mergers. In the hierarchical scenario, 



minor (not equal masses) and major (equal masses) merg- 
ers of low mass galaxies into ones of larger masses is 
the paradigm mechanism by which galaxies acquire di- 
mensions and masses that we observe. There are many 
papers and even text books try i ng to resolve t his is - 
sue, e.g. iBinnev fc Tremaind (|2008l ); iHopkins et all (|2008T) ; 
iNaab et al.l (|2009h to mention a few. The formalism in use 
here is taken from Bernardi (2011, talk delivered at the Uni- 
versity of Padova, private communication). The process can 
be accompanied by star formation (so-called wet merger) 
which would spread the age and the metallicity of the stel- 
lar populations hosted by galaxies. This may cause some 
difficulties with the ETGs whose stars seem to be old with 
a rather narrow age range. We will come back to this later 
on. In alternative, the mergers can occur free of star for- 
mations (in this case they are referred to as dry mergers). 
Wet or dry, by how much a merger between two bodies of 
masses Mi and M2 and radii Ri and R2 can increase the 
mass and size of the daughter galaxy? For a pair of interact- 
ing virialized galaxies made of DM, we may assume that the 
total energy is given by the sum of kinetic and gravitational 
potential 



E — Eyir -(- Eint — Eyir — Kyir -(- Wvir 

where E V i r is the total virial energy of the two galaxies and 
Ei„t is the energy of interaction. This can be neglected with 
respect to the first one so that the total energies before merg- 
ing is 



Ei = 



Mia 2 M 2 a 2 GMj GMj 
_ " : ~Ri ~R~2~ 



2 z jxi ri2 

where the a's are the velocity dispersions. After merging the 
total energy of the system is 



E f = 



(Mi + M 2 )a 2 f G(M 1 +M 2 ) 2 



R f 



In a dissipation-less process, during the merger the total 
energy is conserved, so Ef — Ei. In the case of a major 
merger represented by Mi = M2 = Mi and M/ = 2M^ we 
get 



G{2Mj) 
R t 



2 G(2M f ) 



Rr 



which means that the mass and radius are doubled, whereas 
a remains unchanged. In the MR-plan e the displacement 
vecto r has slope A log R/ A log M — 1 Jvan Dokkum et al.l 
l2010h . 

In the case of a minor merger represented by Mf — (1 + 
f)Mi, where / is smaller than 1, the Virial Theorem (2K = 
—W with obvious meaning of the symbols) yields Mi Ri<r 2 , 
from which we get 

Rf°f = (1 + f)Ri°i = (i + f) 2 ^^/(i + f) 

When / << 1 we obtain 

M f = (1 + f)Mi ~ (1 + 2f)R l a 2 (l - f) 
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i.e. the change in size is larger than the change in mass 
and also a decreases. The displaceme nt vector has slope 
A log R/A log M > 1 likely close to 2 (jvan Dokkum et all 
l2010h . 

There is a final remark worth being made. The above 
formalism has been developed looking at the sole DM com- 
ponent of the interacting galaxies. How the BM component 
will behave during the merger event neither is completely 
understood nor enough simulations are available to make 
solid predictions. Most likely long trails of DM and BM will 
develop, thus likely removing significant amounts of mass of 
the two components. The effect can be particularly relevant 
for BM because of the many implications for subsequent star 
formation (if any). In any case, mergers alone cannot be re- 
sponsible of the well behaved MRR, but would simply add 
some dispersion on it. 

Puffing up and galactic winds. Rapid stripping or ejection 
of BM from galaxies may puff them to larger d imensions. 
The idea goes back to lBiermann fc Shapiro ( 1979) to explain 
the formation of SO galaxies from disc galaxies. Recently 
it has been reconsidered by iRagone-Figueroa fc Granatd 
(|20 1 lT ) to explain the observational evidence that most mas- 
sive ETGs at redshift z > 1 exhibits sizes smaller bay a 
factor of a few than local galaxies of the same type and 
mass. They investigate the effect of BM loss triggered either 
by QSO/starburst/driven galactic winds or quiet gas resti- 
tution and ejection by stars at the end of their evolution, 
both concurring in a galaxy to eject a sizeable fraction of 
the BM mass. The evolution of the galaxy depends on the 
ejection timescale compared to the dynamical timescale. If 
the former is short compared to the latter (fast ejection) the 
ratio of the final to the initial radius is 

R f _ (Mf/Mj) 
Ri 2{M f /M l ) - 1 

and if Mj /Mi < 0.5 the system gets unbound and dissolves. 
If the ejection time scale is long then the conservation of 
adiabatic invariants yields 

Rj_ _ 1 

Rz ~ (Mf/Mi) 

Thus, if the ejection is fast the size increase is more effi- 
cient, if the ejection is slow the system remains bound in- 
dependently of Mf/Mi. In this scheme the compact objects 
should transform in less massive and larger systems. The 
process could occur in two steps: a very fast one at the be- 
ginning of the star forming history a slow one during the re- 
maining life. Indeed our NB-TSPH models partly follow this 
scheme, and partly not. The models have intense episodes 
of star formation and significant galactic winds but on the 
average the general trend is toward larger star masses and 
nearly constant radii. Whether this mechanism can account 
for the general distribution of galaxies on the MR-plane is 
hard to sayQ 

8 The discussion bv lChiosi fc Carrar o (2002) on this issue is still 
valid, and it is strengthened by the new NB-TSPH mod els used 
in this study. In their analysis, IChiosi fc Carrarol (120021) showed 
that the discrepancy between the slopes of the MRR traced by 
models and data for ETGs cannot be ascribed to mergers, because 



Geometrica I e ffects iGuo et all (|2009T ) and 
van Dokkum etafl (|201dh have investigated another 
interesting possibility: the MRR, at least for the most 
massive galaxies, is driven by a systematic variation of 
the Sersic index ng parameterizing t he surface density 
profile s with the redshift. According to Ivan Dokkum et al.l 
in the context of the inside-out growth of a galaxy 
by accumulation of material, the variation of the effective 
radius R e (50% of the light) and R1/2 in turn, the Sersic 
index, and the galaxy mass is 

^ OS f; a 3.561og(n + 3.09) - 1.22 
dlogM 6V ' 

accurate to 0.01 dex for 1 < n < 6. This means that 
the radius increases approximately linearly with the mass 
if the projected density follows an exponential law, but that 
oc M l s for the de Vaucouleurs profile with ns = 4. This im- 
plies that strong evolution in the measured R e is expected 
in all inside-ut growth scenarios irrespective of the physical 
mechanism, unless the density profiles are close to exponen- 
tial. In our models the Sersic index varies from ns — 2 for 
the HDMM case to n s = 4 for the HDHM one. Therefore 
the expected variation of the slope f^^ff due to geomet- 
rical effects goes from about 1 to about 1.8. This implies 
that the displacement vector changes slope according to the 
Sersic index of the underlying galaxy. The slope increases 
from R e oc M to R e oc M 2 . Indeed looking at the mean 
locus of galaxies in the MR-plane, the slope increases mov- 
ing toward the top right corner of the plane, in other words 
the band tends to bend toward higher radii as the mass in- 
creases, passing from 0.5 for the low mass galaxies, to 1 for 
the i ntermediate mass o n es, to about 2 for those at the top 
end. Ivan Dokkum etaTl l|2010t ) made use of the above re- 
lation to explain the steep slope of the MRR for galaxies 
measured at different redshifts (from z = 2 to 2 = 0) whose 
radii and masses seem to increase according to R e oc M 2 : 
0.45 < log ,R e < 1.1 and 11.15 < M(M ) < 11.45. At high 
redshift (z ~ 2) these galaxies have smaller effective radii 
and profiles that are closer to exponential, whereas at low 
redshift (2 ~ 0) the radii get normal and the profiles ap- 
proach the de Vaucouleurs law. It seems that the size and 
hence the density is correlated with the Sersic profile which 
in turn drives the displacement vector in the MR-plane. 
However, it is worth recalling that exactly the same trend 
of the MRR is also shown by the CGS in the mass range 
of interest, thus partially weakening the argument based on 
geometrical effects. Likely both concur to shape the top end 
of the MRR. 



6 THE COMPACT OBJECTS AT HIGH 
REDSHIFT 

As already anticipated the existence of compact, massive 
galaxies is an established fact which seems to be rela- 
tively frequent. Although the detection of compact galax- 
ies at high redshift may be affected by observational 



in this case the models would fall even farther out with respect 
to existing observational data. 
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Figure 12. The logfRjyg] versus log[M 3 ] relation for compact 
galaxies (open ci rcles). The green dots are the normal ETGs of 
the HB sample of Bcrnardi ct al. I J201C|) plotted here for the sake 
of reference. The squares are the twelve reference models together 
with their linear fit (black line), whereas the filled circles are 
the ancillary models with different efficiency of star formation. 
The thin straight lines are the theoretical MRR on which proto- 
galaxies of given DM over-density settle down at the collapse stage 
as a function of t he red shift. The lines are given by eqn. J4} taken 
from iFan et ah and correspond to m = 10 and f a = 1; 

the color-code indicates different redshifts (black: 2 = 0; blue: 
2 = 1; red: z = 2.5; green: 2 = 5; cyan: z = 20). The compact 
galaxies could correspond to models in which star formation has 
been inhibited while the gas was falling into the gravitational 
potential well of DM+BM toward higher and higher densities. 
The models are described in Appendix A and Section 3. 



biases (e. g.. iMancini et al.l 120091 ). extremely deep imag- 
' 20101 ). lensing (e.g., lAuger et alj 
Newton et al.l l201ll ). available measurements of 
very high velocity dis p ersion s (e.g.. | Cenarro fc Truiillol 
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20091: ICappellari et al.1 120091 ; Ivan Dokkum et al.l 120091 

van de Sande et alj|201ll), and m easurements of the surface 



brightness i Szomoru et aljfeoill ). are confirming the i r com - 
pactness. Furthermore, according to ISaracco et all (|201ll ) 
compact galaxies tend to complete their star forming activ- 
ity at z > 5, whereas larger galaxies of similar mass tend to 
have a younger stellar content suggesting that the structure 
of a galaxy may be related to their star for mation activity 
(see lMosleh et aill201ll ; lshankar et al.|[201ll . and references 
therein) . 

The sample of compact galaxies we are going to con- 
sider is gathered from different sources and spans an ample 
range of redsh i fts. D at a for low ma s s obje c ts is taken from 
ICimatti et~ai1 ll2004l), iDaddi et all (120051), iMaraston et all 
(I2006T) . Ffbft et al.l (l2007h IZirm et al.l ll2007l). ICimatti et all 
ihoOSfi. iMcGrath et al.l (l200Sh. Ivan Dokkum et al.l l|2009h . 
GadottiliEoogj), iFisher fc Drorvl l|201ll l. and see also 



respectively. They define a mean relation roughly repre- 
sented by 

logTJi/a = 0.59 log M s - 6.51 
It is worth emphasizing that this rel ation runs nearly parallel 



It is worth emphasizing that this relation runs nearly parall 
to that drawn by the SDSS data of iBernardi et al.1 (2010) 
mas s 
200| ) 



On the h igh 
Bczanson ct al. 
Mancini et al. 



(20 



0) 



mas s ra nge we have data from 
Damianov et all [|2009h . 
Longhetti et al.l ( 20071 ). and re- 



centlv iTiret et al.l ( 20111 ). The estimated masses and radii 
(in the same units) span the intervals 11 < log M s < 12 and 
0.4 < \ogRi/2 < 1-3. They define a mean relation given by 

logi? 1/2 = 1.07 log M s - 11.95 

Remarkably the slope now is nearly twice as much as that in 
the lower mass interval and inte rsects the region of the most 
massive and widest ETGs in the IBernardi et all l|2010h sam- 
ple. For more details on the subject and the current explana- 
tion for both compact and extended objects a t the same red- 
shift ( similar luminosity and hence mass) see lShankar et aD 
(|201lf) . 

The position of the compact galaxies on the MR-plane 
is shown Fig. [12] together with normal ETGs. Even if a for- 
mal comparison between the two samples cannot be made 
because the data are no fully homogeneous, yet an idea of 
their relative position can be inferred. Concerning compact 
galaxies, there are at least three questions to be answered: 
How can these compact (and high mass) objects exist at 
large redshift? How do they evolve? Which are their coun- 
terparts in the local Universe? 

The position of these objects on the MR-plane and com- 
parison with the iso-density lines would imply a rather high 
redshift of formation to which a higher initial density corre- 
sponds. The upper limit is of the order of z ~ 20 — 25. There- 
fore, there are two issues to address: can galaxies with the 
mass currently assigned to compact galaxies exist at these 
redshifts? Is their position indicating the initial density or 
what else? 

Compact galaxies with estimated star mass in the 
ra nge 10 10 to 5 x 10 12 - see for i n stance the data 
bvlGadottil ll2009Ti. iBezanson et aD (I2009TI. Damianov et al. 
(|2009h . IMancini et al.1 ||2010|) . IFisher fc Drorvl l|201ll ). and 
recently ITiret et al.l ( 201 lh - likely belong to systems made 
of DM + BM of significantly bigger mass that can be esti- 
mated from Mom = mM s with m = 15 -j- 25. We adopt 
here m = 1 over the mass r ange 10 10 < M s < 10 13 . 
According to iLukic et al.1 ([20071 ) . the probability of occur- 
rence per (Mpc//i) 3 for galaxies with total mass in the range 
(10 1 " 



Shankar et al.l (|201ll ). The masses (in Mq) and radii (in kpc) 
fall in the ranges 9 < logM s < 11 and 0.3 < logi?i /2 < 1.3, 



10 )M@ goes from a few 10 at £ — 1 to be- 
low 10 at z ~ 16, whereas for objects in the mass range 
(10 12 - 10 13 )M Q the probability falls from a few 10~ 3 at 
z — 1 to below 10 -5 at redshift z ~ 11. These probabilities 
albeit small are still significant for a total observed volume 
of about 10 s Mpc 3 . 

Similar considerations can be made for the objects in 
the lower range of mass even though they have a larger prob- 
ability of coming into existence at high redshifts. So massive 
galaxies with very high initial density (high formation red- 
shift) are possible. 

Looking at the position of the ancillary models on the 
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MR-plane of Fig. 1121 two possible explanations can be put 
forward: i) very high redshift and hence very high mean 
initial density; (ii) more normal redshift and hence initial 
over-densities but intrinsically less efficient star formation. 
This is what indeed is shown by our ancillary models which 
show that in both cases more compact objects are to be ex- 
pected. The case of high redshift, high initial densities is in 
a sense following the general rule. It could correspond to the 
case generating globular clusters and/or bulges of galaxies. 
On the contrary, in the current view of galaxy formation, 
the role played by the star formation threshold density has 
never been fully highlighted and therefore deserves some at- 
tention. If for any reason star formation is inhibited at the 
gas densities that were good for normal galaxies, the pri- 
mordial gas can collapse in the potential well of DM+BM 
to higher densities and smaller radii in turn. When star for- 
mation begins the dimensions of the galaxy get frozen ac- 
cording to our simulations. If the sequence of compact galax- 
ies is generated when a new density threshold is reached by 
the BM inside a DM halo, we expect that it should reflect 
the mean behaviour of the sequence of normal galaxies. We 
have already argued that this latter mirrors the shape of 
the CGS on the MR-plane. Indeed the sequence of compact 
galaxies changes its slope from 0.5 to 1 and even higher as 
the mass goes toward the upper limit. Mergers and mass 
ejection episodes may occur thus inflating the system as al- 
ready reported above and therefore blur the ideal sequence. 
This explanation in terms of the threshold density for star 
formation does not contradict the suggestio n that these ob - 
jects could be the bulges of bigger galaxies (|Graharj20ll . 
Finally, there is another parameter to consider, i.e. f a . Ac- 
cording to eqn. Q, at given total mass, redshift and all re- 
maining parameters, the -R1/2 radius scales with (1.5//o-) 2 . 
Passing from the reference case / CT = 1 to /(r=2, the radius 
goes down by a factor 0.5. Therefore, given the position of a 
compact object on the MR-plane the associated initial red- 
shift would turn out to be about a factor of two smaller 
than before, thus alleviating the initial redshift required to 
explain the position of compact galaxies on the MR-plane. 
The subject is left to future investigation. 



7 SUMMARY AND CONCLUSIONS 

Before drawing some general conclusions about the results 
we have obtained, a few remarks are appropriate. First of all, 
our NB-TSPH galaxy models can be classified as early hi- 
erarchical or quasi monolithic because only during the very 
early stages there exists a hierarchy of substructures merg- 
ing together into a single object. The formation of a mas- 
sive galaxy is com pleted within the first two Gyr or so (see 
iMerlin et ai1l2012l . for more details). Second, the recent dis- 
coveries of high mass galaxies (M s ~ IO^Mq) at high red- 
shift (2 ~ 7 — 8) do not contradict theoretical predictions in 
which the existence of massive galaxies is possible even at 
high redshifts. Indeed, looking at the MR-plane (Fig. [9]), the 
iso-density line corresponding to the collapse redshift 2 = 5 
intersects the CGS at a mass M s ~ 5 x 10 1 Mq, so that 
there is ample room for high mass galaxies. Third, the galaxy 
models in use clearly fol low the so-ca l led do wnsizing pattern, 
originally predicted by ICowfe et all (1 19961), later confirmed 
by the NB-TSPH models of IChiosi fc Carrarol ||2002T ). and 



more rece ntly amply suppor ted by a plenty of observational 
data fsee lBundv et al.ll2006l . for a recent discussion). They 
find that over the probed redshift range, (i) there is a mass 
limit above which star formation appears to be quenched; 
(ii) the mass limit decreases with time by a factor of 3 across 
the redshift range; (iii) the process for quenching star for- 
mation must, primarily, be internally driven. W e point out 
that the arguments given bv lBundv et al. I l|2006h do not con- 
tradict the early hierarchical or quasi monolithic view in 
which galaxies evolved quiescently after the initial period 
of strong mass assembly and star formation. The time in- 
terval span by this early activity was indeed rather short, 
the maximum duration being the time interval elapsed from 
t(zj) to t(z = 1.4). If they later suffered one or more merg- 
ers, most likely these were not accompanied by star forma- 
tion with no consequences on their stellar content. In other 
words, an early-hierarchical scenario cannot be ruled out and 
most likely took place indeed. It is wise to say that early- 
hierarchical, quasi-monolithic and more rarely pure hierar- 
chical schemes all concur to shape the present day popula- 
tion of ETGs; however, the hierarchical scheme cannot be 
the sole dominant channel of galaxy formation. 

In this context, the main results and conclusions of this 
study are: 

• The observational distribution of galaxies on the MR- 
plane does not coincide with locus expected for haloes and 
companion baryonic galaxies on the basis of their mass and 
radii at assigned mean density. Rather, we think that it re- 
sults from two complementary mechanisms: 

(i) As the Universe expands, haloes of higher and higher 
mass collapse. The radius of an halo is proportional to its 
mass according to the relation Rdm oc (iso-density 
line). The proportionality factor depends on the cosmolog- 
ical background and the virialization redshift because sys- 
tems that virialize first have larger mean initial densities 
(the Universe was denser when they collapsed). In other 
words, the initial density of an halo is defined by the back- 
ground density at the moment of its virialization. Inside 
these haloes, the baryonic matter collapses and soon or 
later is turned into stars at a suitable rate, fixing the ratio 
m = Mdm/Ms. So there is a manifold of MRRs described 
by .R1/2 an d M s of a BM galaxy parameterized by the initial 
density, i.e. the redshift of the initial collapse. If the forma- 
tion of a BM galaxy inside a DM halo is nearly dissipation- 
less, the process occurs at constant gravitational potential 
energy per unit mass of the DM and BM components. In 
general, massive galaxies follow this rule, so that the parent 
DM haloes and the daughter BM galaxies lay on nearly par- 
allel MRRs, whose slope is close to 1/3. At decreasing total 
mass of a galaxy, the dissipative processes become more and 
more important so that the MRR slope decreases from 0.333 
for galaxies with total mass 10 13 Mq to about 0.2 for a to- 
tal mass of 10 9 Mq . This latter slope nicely coincides with 
the observational mean value of the MRR for dwarf galaxies 
which are more sensitive to the complex physics of baryons 
(cooling, heating, galactic winds, etc.). 

(ii) The cosmology background fixes the "upper mass" of 
the haloes collapsing at any redshift, i.e. the upper mass 
end of the mass distribution. Assuming a value for the to- 
tal number density of haloes, we derive from the number 
density-plane the CGS to be plotted on the MR-plane in 
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coordinates R1/2 vs M s . The CGS is the locus of the "typ- 
ical masses", function of the redshift, limiting on the high 
mass side the manifold of theoretical MRRs. For 10~ 2 haloes 
per (Mpc//i) 3 the CGS fits astonishingly well the high mass 
boundary of the observational MRR and nicely reproduces 
its slope. Furthermore, the slope of the CGS and that of the 
observational MRR coincide with the relationship predicted 
by the top-hat spherical dissipation- less collapse for primor - 
dial fluctuations originally propo s ed bylGott &: Reesl rtl975Tl 
and fu rther developed by iFaberl (1 19841 ) and Burstein et alj 
This cannot be a mere coincidence. 

• The observational MRR and simulations of the MR- 
plane based on the number densities of haloes n+(MoM, z) 
are in close agreement. The lower boundary of the MRR 
strip is set by the ever increasing typical mass of the under- 
lying mass distribution of the collapsing DM haloes, i.e. the 
CGS. In other words, at any given value of the redshift, i.e. 
along the corresponding MRR, only values below the typ- 
ical mass are permitted. Haloes with higher masses are so 
improbable that the region of the MR-plane below the CGS 
is simply void of galaxies with exception of the so-called 
"compact galaxies", for which a plausible explanation not 
in contradiction with the previous one can be found. The 
upper boundary of the strip is set by the slope of growth 
function of a certain halo mass becoming negative below a 
certain value of the redshift (haloes of that mass start de- 
creasing) plus other physical processes of minor importance. 

• Within the mass interval of existence for haloes formed 
at a given redshift, those with the typical mass form a bary- 
onic galaxy inside closely following the ideal case of the 
dissipation-less collapse. This latter point and the consid- 
eration made above eventually determines the zero point 
and the thickness of the MRR. 

• The low mass end of the simulated MR-plane is over- 
populated by small objects, with respect to the observational 
data. As explained above, this is not a problem, and is in- 
stead consistent both with theoretical expectations and a 
rough counting in the very local Universe (the only region 
of the cosmos in which we can be confident to really obtain 
complete catalogues!). On the other hand, the high mass 
end of the MRR contains galaxies whose total mass falls 
in a risky mass range, more typical of the mass of groups 
or small clusters than of individual objects. Indeed, we do 
not observe galaxies more massive than say M a ~ 1O 12 M 
(the associated halo mass would be ~ 1O 13 M0). In contrast 
we find haloes with masses up M ~ 10 14 - 1O 15 M which 
in turn are hosting groups and clusters of galaxies. In this 
case we may invoke the assumption of the one-to-one cor- 
respondence between halo mass and baryonic galaxy inside, 
and argue that haloes with mass above a given threshold 
value do not contain a single ultra-massive galaxy (with the 
possible exception of the CD galaxies in clusters for which 
recent dry mergers are likely responsible)) but rather host a 
few smaller galaxies. Correcting for this effect would likely 
yield the right number of galaxies also in this region of the 
MR-plane. 

• According to our simulations of the observational MR- 
plane at each redshift the percentage of galaxies undergo- 
ing mergers should be small, otherwise the theoretical MRR 
would not agree with the observational one. More than this 
we cannot say with the present analysis. This is an unex- 
pected important result deserving thorough investigation. 



• Finally, we like to mention that the present analysis 
essentially confirms the suggestion advanced long ago by 
IChiosi fc Carrard (2002) on the base of a much smaller 
sample of observational data and much simpler NB-TSPH 
galaxy models. 
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APPENDIX A 

NB-TSPH MODELS OF ETGS 

iMerlin et af] (|20ld . l2012l ) with aid of the parallel NB-TSPH 
code EvoL produced a number of models for ETGs with 
different mass and/or initial density (twelve cases in total). 
The models are followed from the epoch of their detachment 
from the linear regime, i.e. Zi > 20, to a final epoch (redshift 
Zf) varying from model to model (however with Zf < 1). The 
simulations include radiative cooling down to 10 K (with sets 
a suitable artificial pressure floor to avoid spurious clumping 
of particles), star formation, stellar energy feedback mod- 
eled with a novel method (the Almost-Zero-Mass Particle 
method), re-ionizing photo-heating background, and chem- 
ical enrichment of the interstellar medium. The reader can 
refer to the cited papers for all the details on the method 
and the code in use as well as on model results. 

The assumed cosmology is the standard A— CDM, with 
#0=70.1 km/s/Mpc, flat geometry, fi A =0.721, o- 8 =0.817; 
the baryonic fraction is ~ 0.1656. 

Each NB-TSPH model initially consists of ~ 60, 000 
DM particles plus an equal number of gas particles, a frac- 
tion of which are turned into stars during the simulation. A 
gas particle is eligible to form stars if its density exceeds a 
threshold value (p s / = 5 x 10~ 25 g cm -3 ) and its velocity 
divergence is negative; then, a Monte-Carlo stochastic se- 
lection is adopted to decide which particles effectively turn 
into stars. 

Each model starts as a sphere of DM plus gas particles, 
cut from a wider cosmological simulation in which density 
fluctuations exist and the at the center is a constrained cen- 
tral peak of given density contrast of the cosmological tis- 
sue. The central density peak has a given total mass Mr, 
sum of the DM and BM components, ranging from ~ 10 9 to 
~ 1O 13 M0. This is the proto-halo of our model galaxy inside 
which stars will be formed at later times. The cosmological 
simulation provides the initial positions and velocities of all 
the particles in the prot-halo. A minimal amount of rotation 
and an outward radial component meant to mimic the ex- 
pansion of the Universe are added to all the particles. The 
initial conditions are set in such a way that each model is 
a re-scaled version of a single reference proto-galactic halo, 
with different total mass and/or initial over-density. In this 
way, the attention is focused specifically on the role played 
by different initial masses and densities rather than by local 
inhomogeneities. The properties of all the twelve haloes are 
listed in Tabled 

The proto-galactic haloes are then followed through 
their early stages of expansion following the Hubble flow, 
the turn around, and the collapse. The redshift at which the 
collapse begins varies from model to model and inside the 
same model from the center to the outer regions. In general 
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the collapse occurs during the redshift interval 4 > z > 2, it 
starts first in the central regions and gradually moves out- 
wards. The collapse is complete at redshift z ~ 2. 

All the models develop a central stellar system, with 
a spheroidal shape. They confirm the correlation between 
the initial properties of the proto -haloes and their s t ar for - 
mation histories already found bv lChiosi fc Carrarol <|2002h . 
Massive haloes (Mr — 1O 13 M0) experience a single, in- 
tense burst of star formation (with rates > lO 3 M0/yr) at 
early epochs, consistently with observations, with a less pro- 
nounced dependence on the initial over-density; intermedi- 
ate mass haloes (Mr — 10 11 Mq) have star formation his- 
tories that strongly depend on their initial over-density, i.e. 
from a single peaked to a long lasting period of period of ac- 
tivity with strong fluctuations in the rate; finally small mass 
haloes (Mt — lC^M©) always have fragmented histories, re- 
sulting in multiple stellar populations, due to the so-called 
"galactic breathing" phenomenon. 

The models have morphological, structural and pho- 
tometric properties comparable to real galaxies, in general 
closely matching the observed data; there are minor discrep- 
ancie s that are likely of numerical origin (see I Merlin et al.l 
120121 . for all details). 

These models can be classified as early hierarchical be- 
cause they undergo repeated episodes of mass accretion of 
sub-lumps of matter inside the original density contrast in 
very early epochs and essentially evolve in isolation ever 
since. At the present time, the models closely resemble real 
galaxies. All this leads us to conclude that the evolutionary 
scheme early aggregation followed by isolation (thus much 
resembling the classical monolithic scheme) is able to ex- 
plain many of the observed features of ETGs, particularly 
the complicated and different star formation histories shown 
by haloes of very different mass, without invoking major, 
late mergers (i.e. the classical hierarchical scheme). 

The reference models are calculated adopting a star for- 
mation efficiency e 3 f — 1. This value is larger than current 
estimates fro m observational data in the local Universe, i.e . 
e sf ~ 0.025 jLada fc Ladal 120031 ; iKrumholz fc Tanl l2007t) . 
and theoretical considerations on the global star-to-total 
mass ratio in galaxies suggest e s f < 0.1. However, adopt- 
ing a high value of e s f allows for a strong reduction of the 
computational time, while preserving the bas ic properties 
of th e models (for a complete discussion see I Merlin et al.l 
120121 ). 

For the purposes of this study, using the same numerical 
code and approach for the initial conditions, we also produce 
a group of ancillary models, whose initial parameters and key 
results are listed in Table [3] These models are calculated to 
explore the consequences of much higher initial density con- 
trasts and/or lower star formation efficiencies e s f. The the 
effect of the high er initial density i s alrea dy known from 
th e old calculati o n IChiosi &: Carrard |2002j) and the models 
by iMerlin et~al] (|2012l ). Galaxies of the same mass will be 
shifted on the MR-plane to smaller radii. The effect of star 
formation efficiency is less certain. As expected by chang- 
ing (decreasing) the efficiency e s f star formation is delayed 
or even inhibited. The ga continues to flow into the gravita- 
tional potential well till the threshold density for star forma- 
tion to occur is reached and /or s sufficient number of stars 
are formed, the newly born galaxy has much smaller dimen- 
sion with respect to the corresponding object with higher 



efficiency of star formation. All the ancillary models are cal- 
culated limited to the very early evolutionary stages, to show 
in the MR-plane the initial position of a model galaxy in 
which the gas content has reached densities much higher 
than the formal mean background density fixed by cosmol- 
ogy. The parameter f$ indicates the factor by which the 
initial density of same model in the first group (Table [2]) is 
scaled. Of course by doing this ones has to recompute the ini- 
tial redshift zi and all other parameters using the standard 
procedure based on large scale cosmological simulations. As 
already said in these models we also change the efficiency of 
star formation e s f as indicated in column (3) Table In a 
few models, indicated by e 3 f(z), the efficiency of star forma- 
tion increases with metallicity Z of the gas content accord- 
ing to e sf = MIN(1., io.(- 5 '°9( z )+ 1 )) on the notion that a 
metal-rich gas finds it easier to collapse and form stars. The 
efficiency goes from e 3 f =0.1 for Z=0.0001 to e s f = 1 for 
Z=0.01 (close to the solar value). No special meaning must 
be given to this relation, it is simply meant to evaluate the 
effect of an efficiency of star formation increasing with the 
metallicity. In any case, this effect plays a marginal role on 
the position of the model galaxies on the MR-plane, see the 
entries of Table [3] The effects of higher initial density and 
efficiency of star formation passing from e s f = 1 to 0.1 or so 
are of paramount importance and cannot be ignored. 
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Table 2. The twelve reference models. Left to right: total initial mass My = Mbm + MgM [ m units of 1O 12 M0], corresponding (initial) 
gas mass M g [in units of 10 12 AfQ], initial redshift z%, mean halo over-densities [8p — l]j at the initial rcdshift, initial proper physical 
radius of the halo [in kpc] , redshift of the last computed model z t , corresponding age t t [in Gyr] , virial radius of the whole system R v i r 
[in kpc], half-mass radius R\/2 at Zf, and the total stellar mass M 3 at final redshift [in units of 10 12 Mq]. 
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Table 3. The ancillary models. The meaning of the symbols is as follows: Model is the two-letter string identifying the model according 
to the mass: MM for intermediate mass galaxy 2.69 X lO n M and LM for the low mass case 4.17 X 10 9 Mq; f$ is the multiplicative factor 
of the initial over-density, in other words the starting over-density of the simulation is a factor fg higher that the standard over-density 
currently assumed for the reference model of the same mass; e 3 f is the dimensionless efficiency of the star formation rate, the symbol 
e 3 f(Z) means that the efficiency is supposed to increase from e a f = 0.1 for Z=0.0001 to e s f = 1 for Z=0.01 (close to the solar value). 
All other symbols have the same meaning as in Table [2] 
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181 


0.18 


6.3 


49.0 


0.05 


2.6 


0.79 


0.020 


MM 


20 


1 


0.269 


0.0445 


181 


0.18 


6.3 


30.4 


0.10 


4.9 


1.26 


0.054 


MM 


20 


0.1 


0.269 


0.0445 


181 


0.18 


6.3 


56.0 


0.04 


2.1 


0.08 


0.0050 


MM 


20 


0.1 


0.269 


0.0445 


181 


0.18 


6.3 


49.0 


0.05 


2.6 


0.08 


0.0074 


MM 


20 


e e .f(Z) 


0.269 


0.0445 


181 


0.18 


6.3 


59.0 


0.04 


2.0 


0.07 


0.0036 


MM 


20 


e sf (Z) 


0.269 


0.0445 


181 


0.18 


6.3 


46.0 


0.05 


2.9 


0.09 


0.0083 


MM 


15 


1 


0.269 


0.0445 


140 


0.18 


8.2 


30.4 


0.10 


4.5 


1.08 


0.026 


MM 


12 


1 


0.269 


0.0445 


117 


0.17 


9.8 


30.4 


0.10 


4.1 


0.92 


0.011 


MM 


12 


1 


0.269 


0.0445 


117 


0.17 


9.8 


18.8 


0.20 


7.4 


0.98 


0.015 


MM 


12 


e sf (Z) 


0.269 


0.0445 


117 


0.17 


9.8 


36.0 


0.08 


3.1 


0.07 


0.0032 


MM 


12 


e sf (Z) 


0.269 


0.0445 


117 


0.17 


9.8 


30.4 


0.10 


4.2 


0.10 


0.00575 


MM 


5 


1 


0.269 


0.0445 


68 


0.14 


16.7 


7.0 


0.70 


19.2 


5.75 


0.019 


MM 


5 


1 


0.269 


0.0445 


68 


0.14 


16.7 


4.1 


1.50 


34.3 


6.76 


0.057 


LM 


20 


1 


0.00417 


0.000691 


112 


0.16 


2.6 


25.3 


0.13 


1.2 


0.32 


0.00016 


LM 


20 


1 


0.00417 


0.000691 


112 


0.16 


2.6 


18.8 


0.20 


1.8 


0.31 


0.00017 


LM 


20 


0.1 


0.00417 


0.000691 


112 


0.16 


2.6 


25.4 


0.13 


1.3 


0.02 


0.000091 



